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ABSTRACT 

The ability of a CFD technique to evaluate the performance of a heat 
shield separation system has been demonstrated. The evolution of 
the differential pressure across the heat shield surface has been 
estimated. The results of the analysis have been compared with the 
data from the separation test of the heat shield of an Indian launch 
vehicle. The pressure on the heat shield surface was estimated by 
solving the Navier Stokes equations using the Eulerian Lagrangian 
formulation. The results, which show good agreement with the 
observed results of a heat shield separation test on the ground, were 
found to be adequate for design purposes. 


NOMENCLATURE 

D velocity divergence 

f r body force component in radial direction 

f x body force component in axial direction 

/ internal energy 

p pressure 

R radius of heat shield fairing 

t time 

u velocity component in axial direction 

v velocity component in radial direction 

x, y coordinate directions 

p density 

A., |L coefficients of viscosity 

n„, n vv , n, v n 9 stress deviators 


1.0 INTRODUCTION 

The payload fairing in a rocket protects the satellite from aerodynamic 
loading, heating, acoustic vibrations and other adverse environmental 
conditions during launch and the ascent flight. The fairings are 
separated as soon as the vehicle gets outside the atmosphere. 
Convincing demonstration of jettisoning of the heat shield fairing — 
clean separation of the fairing halves, devoid of any possibility of re¬ 
contact with the vehicle or parts of the payload — is an imperative step 
in the development of any launch vehicle. During flight, the separation 


event occurs in near vacuum condition. The stringent requirements on 
the dynamic clearance of modem day fairings coupled with the need to 
keep the weight low, call for design of highly flexible fairing halves 
which undergo considerable deflection during separation. While this 
flexibility and consequent large deflection helps in attainment of large 
dynamic clearances during separation, it also makes the structural 
behaviour of the fairing half to become highly sensitive to minor 
changes in loading on the heat shield surface as, for example, may 
happen if the test were conducted in air instead of vacuum. Though the 
conduct of such a separation test on full scale model necessitates the 
availability of a large vacuum chamber and considerable expenditure, 
the importance of performing a reliable demonstration of the event has 
prompted the designers of rocket vehicles to persist in performing the 
test under vacuum conditions. For example the separation test of the 
flight standard version of Ariane 4 payload fairing (4m diameter and 
9-6m height) was conducted in the Dynamic Test Chamber at ESTEC, 
Netherlands, stated to be the largest vacuum chamber in Europe. A 
detailed description of this test is available in Ref. 1. Similarly the 
separation tests on the payload fairings of the Atlas launch vehicle® 
were conducted in the world’s largest vacuum chamber, at the Space 
Power Facility, Plum Brook, after simulating the altitude of jettison 
event. 

In contrast with the usual method of qualification employed in the 
case of Ariane and Atlas heat-shields, the test on the heat shield of the 
Indian Polar Satellite Launch Vehicle (PSLV) was performed in air and 
the effects of the air were accounted for by analysis. The present study 
is an excellent example of the application of CFD capabilities for the 
solution of a practical problem of interest. A similar approach, though 
complete details are not available, also appears to have been adopted 
by the designers of the Japanese H-II rocket 13 ). While the report does 
mention the measurement of the pressure inside the heat shield 
compartment, it is not clear whether any attempt was made to predict 
the heat shield response on the basis of analysis only or whether the 
experimental pressure data were used as the sole inputs for verifying 
the response. Present results have generated a sufficient degree of 
confidence to enable the designer to extrapolate the results of experiments 
in air to vacuum conditions. Heat shield separation system of PSLV 
was designed without conducting any test in a vacuum chamber, leading 
to large cost savings. Successful separation of heat shield fairings of all 
PSLV missions demonstrated the correctness of the approach. 


Paper No. 2427. Manuscript received 16 June 1998, revised version received 5 June 1999, accepted 13 July 1999. 




436 


The Aeronautical Journal 


September 1999 


2.0 SEPARATION TEST 


2.1 The experimental set up 

Figure 1 shows a sketch of the PSLV heat shield separation test set 
up. The test article consists of a flight quality PSLV heat shield 
mounted on a simulated lower stage interface. The heat shield of 
aluminum isogrid construction consists of a conical upper section 
with spherical nose, a cylindrical middle section and a lower boat- 
tail section which mates with the top of the third stage of the launch 
vehicle. The shield is made of two half fairings joined along the 
longitudinal seam by the help of piston cylinder assembly with rubber 
bellows between the two. The cone-cylinder-boat-tail assembly is 
attached to the lower stage through a mermen band clamp. Up to the 
time of deployment, the heat shield fairing halves are held together 
— the piston being held in place within the cylinder — by shear rivets. 
The explosive cords are inside the bellows, kept within the cavity of 
the assembly. The bellows are bonded to the cylinder to prevent their 
detachment during the separation process. Two bellow systems are 
placed at the two diametrically opposite edges of the fairings. These 
are invariably interconnected at the top to provide redundancy. The 
jettisoning of the heat shield is achieved in two steps. The first consists 
in activating a horizontal separation system which detaches it from 
the parent vehicle. This is followed by lateral separation system 
which throws the fairing halves apart. A single command 
simultaneously activates the pyrotechnic bolt cutter releasing the 
mermen band and the ignition of the explosive in the zip cord. Ignition 
of explosive in the zip chord causes the bellow pressure to increase 


enormously, resulting in the shearing of the rivets and the impulsive 
ejection of the two halves of the fairing. 

Two tests were conducted. The first one took place before the 
effects of the air on separation were clearly understood and the second 
one was repeated with additional instrumentation consisting of: 

1. Pressure transducers to record the pressure inside the heat shield 
compartment and on the surface. 

2. Strain gauges mounted on the surface. 

3. Accelerometers on the heat-shield surface. 

4. Measurement of bellow pressure to compute delivered impulse. 

5. Optical sensors to determine the time when bellow is losing contact 
with the piston. 

The locations of the pressure transducers (cp t to cp l2 ) are shown 
in Fig. 1. The achievement of minimum separation distance was 
confirmed by the use of snap gauges. A detailed description of the 
experiments conducted and the data generated has been presented in 
Ref. 4. 


3.0 ANALYSIS 

In this paper we have evaluated the transient behaviour of the 
compartmental pressure by solving the two-dimensional Navier 
Stokes equations. A Simplified Arbitrary Lagrangian Eulerian 
(SALE) algorithm developed by Amsden, Hirt and Cook< 5 ' 6 > has 
been used to solve the equations after necessary modifications. The 
most important feature of this algorithm is that it uses both Eulerian 
and Lagrangian procedures to solve the governing equations and is 
applicable to flows at all speeds. The procedure uses a finite difference 



are located in the 
outer surface 


Figure 1. Set up for separation test of hield shield. 
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mesh with the vertices that may move with the fluid (Lagrangian), be 
held fixed (Eulerian) or be moved in any other prescribed way. The 
flow variables, pressure (p ), density (p) and internal energy (/) are all 
assigned to cell centres, while the coordinates (x, y) and velocity 
components («, v) are assigned to the cell vertices. The calculations 
necessary to advance the solution for one time step are performed in 
three distinct phases. In the first phase, explicit Lagrangian calculations 
are made, where convective fluxes are zero. In the second phase, 
which is optional, implicit calculations are made to obtain time 
advanced pressure forces. The main advantage of the second phase 
is that it eliminates the numerical stability condition that limits a 
sound wave to travel no further than one cell per time step. The cell 
vertices are moved to their new Lagrangian position after this phase. 
Finally in the third phase, which is also optional, the mesh can be 
relocated in a new position. In this re-zone phase, convective fluxes 
are computed to account for the movement of the fluid between cells 
as the mesh moves. Further details of the SALE algorithm are 
available in Ref. 6. In the present work, we have used the explicit 
Eulerian calculation only. 


3.1 Governing equations 

The Navier Stokes Equations in the axisymmetric form can be written 
as 

Mass equation: 

3p + J_3rpi£ + apv =0 (1) 

3 1 r dx dy 
Axial momentum equation: 


3p« | 1 drpu 2 | 3puv _ dp | 1 3rll M | 
3 1 r dx dy dy r dx 


i!k + rk +p/ . 


dy r 

Radial momentum equation: 

3pv 1 3rpwv 3pv 2 dp 1 dr\{ 311™ r 

Energy equation: 

3p/ 1 drplu dplv „ „ du „ du nFL 

„ 3v „ 3v 
+ H ” 3x + 0 vv dy 

where D is the velocity divergence: 

D _ 1 dru + dv 
r dx dy 


- ( 2 ) 


• • (3) 


(4) 


.(5) 


Velocity components («, v) are either in the Cartesian coordinates (x, 
y) or in the cylindrical coordinate (r, z). When Cartesian coordinates 
are desired, all radii r which appear in these equations are set to unity. 
The fluid pressure (p ) is determined from the equation of state 


P = P{ P’ 7 ) 


• • • (6) 


The stress deviators are defined according to 

n„=2 pg+w) 


3v 

'dy 


n yy = 2p— + XD 


n 0 = CFL^2p- + Xi)j 

(du 3v 


... (7) 


'’■’-'‘UTS 


where X and |l are the coefficients of viscosity connected by the 

Stokes law, viz 

3X +2(1 = 0 ...(8) 


The coefficient CYL is zero for Cartesian coordinates and unity for 
cylindrical coordinates. 

3.2 Idealisation of the problem 

The firing of the explosive cord generates enormous lateral force 
projecting the fairing halves away from each other with a very large 
acceleration. Figure 2 depicts a typical recording of accelerometer 
data obtained during the test on the heat shield. This sudden motion 
in the presence of air gives rise to flow phenomena akin to the flow 
on the two sides of the diaphragm in a shock tube. The pressure 
differential between the two sides of the heat shield plays the primary 
role in diminishing the velocity gained and also in modifying the 
structural response of the heat shield to the impulsive forces acting 
on its edges. During this period the flow field is dominated by the 
waves generated by the sudden motion of the heat shield surface and 
their reflections. 

For the purpose of estimating the compartmental pressure, the 
heat shield configuration (Fig. 1) is idealised as consisting of two 
concentric, infinitely long circular cylinders separating across a 
diametrical plane. We are thus ignoring the effect of the closed nose 
cone, the edge effects at the top and bottom extremities and also the 
presence of a number of components inside the heat shield which 
will modify the reflection of the waves. The radius of the inner cylinder 
corresponds to the radius of the spacecraft and the radius of the outer 
cylinder corresponds to the radius of the heat shield fairing. It is 
intuitively recognised that the pressure differential observed in the 
experiment will be higher than our estimates at the top and lower 
than the estimate at the bottom. The resulting occurrence of a nose 
closing tendency in the early phase of heat shield separation was 
indeed noted in the experiment. 

By taking advantage of the symmetry condition, the computation 
is performed only for the upper half of the semi-circle. 

3.3 Modelling of zipcord acceleration and atmospheric 
drag 

During the first few milliseconds, when the zipcord system is in contact 
with the heat shield, the acceleration of the fairing is mainly dominated 
by the expansion of the bellow. The acceleration is found to reach its 
peak almost instantaneously with the initiation of the separation 
command and fall equally rapidly to zero when the contact with the 
bellow is lost. In the following period, the fairing is subjected to 
atmospheric drag due to the pressure differential between the outer 



Figure 2. Typical acceleration profile. 
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surface and the inner surface. A linear model is employed in the 
simulation for the zipcord acceleration, and the atmospheric drag has 
been calculated at every time step from the computed pressure 
differential between the outer and inner surface of the heat shield 
fairing. 

3.4 Computational grid 

The grid network for the computational domain is shown in Fig. 3. 
The bottom boundary (AB) corresponds to the surface of the inner 
body (payload) inside the heat shield compartment and it is aligned 
with one mesh line. The heat shield surface is also aligned with a 
mesh line (EE) along the radial direction. The outer computational 
boundary (CD) is placed at a distance of five heat shield diameters 
beyond the heat shield fairing surface. The left boundary (AD) is 
placed in the plane of separation and in the plane of symmetry the 
center line boundary (BC) is placed. 

A 38 x 40 grid is taken for the computation. In order to capture 
the complexity of the flow field, the mesh is relatively fine near the 
line of separation and the grid is exponentially stretched in the 
circumferential direction. In the radial direction, a uniform grid is 
employed in the gap between inner body surface and the heat shield 
fairing and in the region outside the heat shield fairing the grid is 
stretched exponentially. 

3.5 Boundary conditions 

Four distinct types of boundary conditions were used in the present 
analysis: 

1. Inflow boundary 

2. Outflow boundary 

3. Symmetry boundary 

4. Solid wall boundary 

In the inflow boundary (AD), not only the velocity components 
but also the density and energy have to be prescribed. A typical 
treatment of the outflow boundary (DC) is either to specify the 
outflow velocities or set the velocity components equal to those 
located one vertex in. This generally works properly for high speed 



Figure 3. Computational grid. 


flow, but for low speed or incompressible flow, the prescription 
requires adjustment if it affects the upstream flows. In the present 
analysis we have used a specified pressure boundary. It uses fictitious 
cells for pressure definition. For a symmetry boundary condition 
(BC), the normal wall velocity must be kept zero throughout the 
calculation. If such a boundary is parallel to the x or y axis, the axial 
velocity is set to zero on a left or right boundary, or the radial 
component of velocity is set to zero on a bottom or top boundary. If 
the wall is slanted or curved, both the velocity components must be 
adjusted to make the flow tangential to the local slope. For the solid 
wall boundary (AB), the velocity components are set to zero regardless 
of wall orientation or curvature. 


4.0 RESULTS AND DISCUSSIONS 

The validation of the software was done by comparing the computed 
results with the experimental results of the PSLV heat shield separation 
test on air. In the beginning a small gap of 5mm between the heat 
shield fairing and the plane of separation was provided to allow the 
flow inside and outside the compartment to interact. This initial gap 
is kept fixed until 5ms when the bellow loses contact with the heat 
shield fairing as observed in the experiment. Afterwards the new 
position of the fairing is calculated after every millisecond from the 
computed velocity and the grid is redefined for that new position of 
the heat shield fairing. 

In order to estimate the variation of pressure drop at various 
locations inside the heat shield compartment, detailed measurements^ 
were conducted. Figure 1 gives the location of the pressure sensors 
for the experiment. Two sensors (cp 2 , cp 5 ) were provided for the outer 
surface of the heat shield fairing. Four sensors (cp w , cp 9 , cp 4 , cp 2 ) to 
determine the longitudinal variation and another four sensors (cp-,, 
cp 6 , cp 4 , cp g ) were provided inside the compartment to determine the 
circumferential variations of temporal pressure behaviour. Two other 
sensors (cp n , cp i2 ) were placed at the R h and K U location at the centre 
of gravity (CG) plane inside the heat shield compartment. 

Figure 4 presents the comparison of measured and computed values 
of the transient delta pressure (kPa) (differential pressure from the 
atmospheric value) at the inner surface at the CG plane. Good overall 
agreement is observed between the two, except for the time of 
occurrence of the pressure recovery. The computed time of occurrence 
of the pressure recovery is about 10ms later than that of the 
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Figure 4. Pressure profile at inner surface at the CG plane. 
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experimental value. The initial overshoot present in the measured 
data during 10-20ms is due to experimental noise. As the command 
for ignition of explosive cord is given only at 19ms, no changes in 
the pressure are expected before that instant. The difference in the 
magnitude of the peak pressure drop between the measured and the 
predicted value is about 20%. 

Figure 5 gives the transient pressure profile at the plane of symmetry 
on the outer surface of the heat shield fairing. Pressure at the outer 
surface keeps on increasing up to 5ms because of zip cord acceleration. 
After 5ms, when the fairing loses contact with the bellow, the pressure 
at the surface slowly comes down to atmospheric level. The small 
oscillation in the pressure in the later part is due to the reflection of 
the waves from the outer computational boundary and also due to 
numerical oscillation. Comparison with the experimental data shows 
a good overall agreement. 

Figure 6 depicts the measured internal surface pressure drop at 
various circumferential locations. It is noted that at any particular 
instant of time, the circumferential variation of the surface pressure 
drop is almost negligible. This may be compared with the computed 
values of the circumferential variation of the pressure drop at the 
instant of 29-39 and 34-31 milliseconds given in Fig. 7. 

Figure 8 shows the comparison of experimental and computed 
pressure drop at the inner surface and the R h location at the plane of 
symmetry. As expected on the inner surface the pressure drop starts 
from 19ms onwards, whereas at the R h location the pressure remains 
constant at the atmospheric level until the pressure waves reach the 
centrebody surface. The time matches well with the time the sound 
wave take to travel the distance of 0-8m (the distance between the 
centrebody surface and the heat shield fairing). This provides a 
check on the correctness of the calculation. The peak pressure drop 
and the overshoot for R h location is slightly more than that of the 
inner surface location. In the remaining part, the behaviour of the 
two curves is almost similar. The comparisons of the experimental 
data for the two locations also show a similar trend. The longitudinal 
variation of the internal pressure drop observed in the test is given in 
Fig. 9. As expected, the pressure drop at the top location (near the 
nose cap) is more than that at the bottom location. It may be recalled 
that the present computation cannot predict this longitudinal variation 
because the formulation assumes two-dimensionality. 

Figure 10 gives the computed acceleration profile of the fairing. 
Up to 5ms, when the fairing is in contact with the zip cord bellow. 




Figure 6. Circumferential variation of internal surface pressure drop 
(experimental). 
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Figure 7. Circumferential variation of internal surface pressure drop 
(computational). 


the acceleration is mainly dominated by the zip cord acceleration. 
When the fairing loses contact with the bellow, it is subjected only to 
atmospheric drag which is caused by the pressure differential between 
inner and outer surfaces. Figure 11 presents the pressure profile at 
the plane of separation. The unsteady nature of the flow field appears 
truly physical, related to the vortex shedding from the edge of the 
fairing. As our interest is primarily related to the surface pressure on 
the heat shield, no attempt was made to capture the vortex motion in 
the computation. This unsteadiness is confined upto a few cells near 
the plane of separation. This is not found to influence the flow field 
in the plane of symmetry. 

In order to analyse the effect of grid size, simulations were also 
carried out for another grid of size 50 x 46. Figure 12 compares the 


Figure 5. Pressure profile at outer surface. 







440 


The Aeronautical Journal 


September 1999 



Figure 8. Radial variation of pressure drop. 



Figure 9. Longitudinal variation of internal surface pressure drop. 

compartmental pressure drop at the R h location for the two grids at 
the plane of symmetry, which clearly indicates the grid independence 
of the results. 

In order to find out the effect of the initial gap between the plane 
of separation and the heat shield fairing in the solutions, computations 
were performed by taking the initial gap of 5mm, 15mm and 30mm. 
Transient pressure profiles at the R h location with different initial 
gaps are shown in Fig. 13. It has been found that with the increase in 
the gap the maximum pressure drop decreases. Moreover the 
overshoot, which is present for the gap of 5mm, is not noticed for the 
other two cases. Based on this observation, we can infer that for a 
proper simulation it is necessary to take the initial gap as small as 
possible. 

For the stability of the scheme the following conditions on time 
step as suggested by Amsden et al ,6 > is observed throughout the 



Figure 10. Acceleration profile. 
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Figure 11. Pressure distribution on outer surface at the plane of separation, 
computation: 

• Courant condition on sound signal propagation is not violated. 

• Fluid cannot move more than approximately one cell width per 
time step. 

• For viscous effect, in every cell momentum must diffuse less than 
one cell width per time step. 

• Time step should not increase more than 5% from the time step 
of the previous cycle. 

• Time step should not be greater than the maximum time step for 
calculation (prescribed). 

SALE algorithm is generally considered to be first order accurate' 6 ). 
A spectrum of computer runs with different mesh, time steps, artificial 
viscosity, convergence criteria, etc, were made to determine if a 
calculated effect is physical or simply a numerical artifice. Apart 
from this type of brute force approach, it has also been checked that 
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Figure 12. Pressure profile at % location (grid independence result). 



the solution does not exhibit large variations over distances comparable 
to cell width. 

The simulations were carried out in the workstation based on Intel 
860 chips. For a moderate grid of 38 x 40, the software take about 
0-18 seconds per iteration. The total number of iterations required for 
calculating the flow field for 100ms is about 60,000. 


5.0 CONCLUSIONS 

To explain the phenomena of the compartmental pressure drop and 
to simulate the aerodynamics of the heat shield separation on the 
ground, software has been developed to solve the two-dimensional 
Navier-Stokes equations by using the Simplified Arbitrary Lagrangian 
Eulerian (SALE) numerical algorithm. The developed software was 


validated by comparing the computed results with the results of the 
full scale heat shield qualification test of the Indian Polar Satellite 
Launch Vehicle on the ground. 

The prediction captures all essential features of the temporal 
behaviour of internal pressure including the time of occurrence of 
peak, overshoot, etc. The difference of the computed magnitude of 
the peak drop to that of experiment is about 20%. A good match has 
been observed between experiment and computation for the pressure 
profile at the outer surface of the heat shield fairing. The computed 
results indicated that the pressure at the R h location was not very 
different from that on the internal surface as observed in the 
experiment. No appreciable variation of internal pressure drop along 
the circumference has been observed in both experiment and 
computation. Understanding of aerodynamic process involved in the 
separation test in the air was found to be very helpful in extrapolating 
the test results to the flight conditions and a successful design of a 
heat shield separation system was made without a test in vacuum 
chamber, leading to a large cost savings. 

It is necessary to extend the program for a three-dimensional case 
to capture all the behaviour of the internal pressure drop in the 
longitudinal direction inside the heat shield compartment. It may be 
noted that the edge effects will introduce a nose closing tendency 
during separation. The presence of such a tendency was confirmed 
by the test data. Time accuracy has also to be introduced to make it 
faster and more efficient numerically. 
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